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Abstract. The southwestern Balkans is inhabited by three endemic water frog species, Pelophylax shqipericus, P. epeiroticus 
and P. kurtmuelleri. Although these species are genetically and bioacoustically distinct from each other, their morphologi- 
cal differentiation is less pronounced as is generally observed in the genus Pelophylax. In this study, we analyzed morpho- 
logical data of 215 individuals of these Balkan endemics in order to evaluate their morphological variability and facilitate 
their identification in the field. Moreover, since P. kurtmuelleri is a disputable taxon, we compared it with closely related 
P. ridibundus from Central Europe. The most pronounced morphological differences expressed by the analysis of variance 
(ANOVA) and the discriminant analysis of the principal components (DAPC) were found between P. kurtmuelleri and 
the two other Balkan species. Pelophylax kurtmuelleri and P. ridibundus, as well as P. shqipericus and P. epeiroticus were 
only slightly differentiated by morphology, which demonstrates that the morphological characters measured do not reflect 
the genetic divergence of the species and are weak indicators of their taxonomic status. Morphological characters suitable 
for species identification include a combination of indices L/CINT, T/CINT, L/T, and DP/CINT (L - snout-vent length, 
CINT-length of the metatarsal tubercle, T - length of the tibia, DP - length of the first toe), and qualitative traits like the 
pattern and shape of the lateral spots, presence/absence of yellow pigment in the flanks and thighs, foot webbing coloura- 
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tion, size and shape of the metatarsal tubercle, and colouration of the male vocal sacs. 
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Introduction 


European water frogs (genus Pelophylax) are a widespread 
group of amphibians and have been intensively studied 
since CARL VON Linné described the first representative, 
Pelophylax esculentus (LINNAEUS, 1758). Modern herpeto- 
logical research has since led to the discovery of up to 22 
species and three hybrid asexual forms, which are distrib- 
uted throughout the temperate and subtropical zones of 
the whole Palaearctic (PLOTNER 2005, FrosT et al. 2020). 
The southwestern Balkans is inhabited by three water 
frog taxa, all of them endemic to the peninsula. The Alba- 
nian water frog, P. shqipericus (Hotz, UzzELL, GUNTHER, 
TUNNER & HEPPICH, 1987), described from the Lake Ska- 
dar in Montenegro, is distributed throughout the lowlands 


of the Adriatic coast of southern Montenegro and western 
Albania (GUNTHER 2004, PLOTNER 2005, VUCIĆ et al. 2018). 
Mitochondrial phylogenetic analyses has shown that this 
species is closely related to P. lessonae (CAMERANO, 1882) 
and P bergeri (GUNTHER, 1986) (RAGGHIANTI et al. 2004, 
AKIN et al. 2010, PLOTNER et al. 2012, DUFRESNES et al. 2017), 
forming the so-called P. lessonae group (cf. PLOTNER 2005). 
The Epirus water frog, P. epeiroticus (SCHNEIDER, SOFIANI- 
DOU & KYRIAKOPOULOU-SKLAVOUNOU, 1984) was described 
from the Lake Ioannina (Greece) and is currently known 
from southern Albania and western Greece (SCHNEIDER 
et al. 1984, SOFIANIDOU & SCHNEIDER 1989, SCHNEIDER & 
HAXHIU 1992, ORUÇI 2010). Its phylogenetic position is not 
clearly resolved but the species is likely more related to P ridi- 
bundus (PALLAS, 1771) and other species from the so-called 
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P ridibundus group (named according to AKIN et al. 2010), 
than to the P lessonae group (cf. LYMBERAKIS et al. 2007, 
AKIN et al. 2010). Another endemic taxon from the Balkans, 
P. kurtmuelleri (Gaypa, 1940) was described as early as the 
19408, but until recently was not recognized as distinct from 
P ridibundus. Research in the early 1990s, however, revealed 
that populations from the southwestern Balkans differed 
from other European P. ridibundus populations on the ba- 
sis of bioacoustics (SCHNEIDER & SINSCH 1992, SCHNEIDER 
et al. 1993) and allozymes (SINSCH & EBLENKAMP 1994, SO- 
FIANIDOU et al. 1994). Currently, the taxon is recognized as 
a separate species (PLOTNER et al. 2012, Frost et al. 2020) 
from, or a divergent evolutionary lineage of P. ridibundus 
(LYMBERAKIS et al. 2007, SPEYBROECK et al. 2020). 
Pelophylax kurtmuelleri lives sympatrically and often 
also syntopically, with P. shqipericus or P. epeiroticus and 
occasionally hybridizes with them (Hotz & UzzELL 1983, 
SCHNEIDER et al. 1984, Horz et al. 1985, SOFIANIDOU & 
SCHNEIDER 1989, Vuci¢ et al. 2018). Ranges of the last two 
species do not overlap and are separated by a 60 km long 
mountain range in southern Albania inhabited only by 
P. kurtmuelleri (PLOTNER 2005, SPEYBROECK et al. 2016). 
However, hybridization between P kurtmuelleri and its two 
endemic counterparts seems uncommon and hybrids are 
rare and likely do not exhibit hybridogenetic reproduction 
as opposed to other regions of Europe (Hotz & UZZELL 
1983, Hotz et al. 1985, SOFIANIDOU & SCHNEIDER 1989). 
Water frogs have relatively uniform morphology, and, 
therefore, their systematics have often relied on genetic and 
bioacoustic characters. Morphology is, however, necessary 
for taxonomists who use morphological characters when 
describing new taxa, for evolutionary biologists studying 
adaptations, phenotypic plasticity, or life-history traits, 


and finally, for field biologists who need to identify species 
directly in the field. Morphological examination of taxa 
currently recognized as P. shqipericus, P. epeiroticus, and 
P. kurtmuelleri was carried out for the first time by Hotz & 
UzzELL (1982), later by SCHNEIDER et al. (1984, 1993) and 
GAVRILOVIĆ et al. (1999). These authors found significant 
differences between the species in several morphometric 
indices and qualitative traits, such as colouration, head 
size, and size and shape of the metatarsal tubercle. 

In this study, we aim to build on previous studies and re- 
vise the morphological characters and morphometric indi- 
ces for three Balkan and one Central-European water frog 
species. Here we benefit from proper identification of each 
individual using genetic markers and specifically 1) exam- 
ine external morphology of P. shqipericus, P. epeiroticus and 
P. kurtmuelleri and determine the quantitative and quali- 
tative morphological traits suitable for their identification, 
and 2) compare morphological variation of P. kurtmuelleri 
with its closely related species, P. ridibundus. 


Material and methods 
Sampling and data collection 


We collected 215 water frogs of the genus Pelophylax from 
13 localities in Albania, Greece, and Montenegro between 
2017 and 2019 (Tab. 1, Fig. 1). Selected localities covered 
the range of the species P. epeiroticus, P. kurtmuelleri and 
P. shqipericus. Since P. kurtmuelleri is closely related to 
Central-European P. ridibundus phylogenetically, we also 
included 89 individuals of this species for comparison. All 
P. ridibundus individuals were sampled during the years 
1997 to 2009 in five localities in western Slovakia. 


Table 1. A list of localities and water frog species (genus Pelophylax) included in the study. 


Coordinates P epeiroticus P. kurtmuelleri P. shqipericus P ridibundus 
Locality Country N E Males Females Males Females Males Females Males Females 
Berat Albania 40.7043 19.9480 - - 1 - - - - - 
Dushnik Albania 40.7594 19.9891 - - - 2 - - - - 
Qazim Pali Albania 40.0497 19.8420 - - 15 9 - - - - 
Ksamil Albania 39.7454 20.0207 - 1 1 1 - - - - 
Mbrostar Albania 40.7743 19.5897 - - 2 3 - - - - 
Nishaj Albania 41.6852 19.5877 - - - - 27 19 - - 
Divjaké - Karavasté Albania 40.9901 19.4975 - - 10 11 7 4 - - 
Poçem Albania 40.5051 19.7036 - - 22 10 - - - - 
Xare Albania 39.7330 20.0529 1 - - 3 - - - - 
Zvërnec Albania 40.5179 19.4007 - - 1 - 3 - - - 
Igoumenitsa Greece 39.5366 20.2038 15 5 4 4 - - = 7 
Ioannina Greece 39.6886 20.8584 9 20 3 1 - - - - 
Virpazar Montenegro 42.2450 19.0915 - = = - - 2 - 7 
Brodské Slovakia 48.6951 17.0089 - - - - - 3 1 
Čunovo Slovakia 48.0331 17.2006 - - - - - - 3 1 
Devín Slovakia 48.1749 16.9766 - - - - - - 15 9 
Šúr Slovakia 48.2268 17.2054 - - - - - - 23 13 
Šaštín-Stráže Slovakia 48.6379 17.1400 - - - - - - 13 8 
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Morphological differentiation of Pelophylax from the southwestern Balkans 


Each individual was measured using a digital calliper 
for five morphometric characteristics: snout-vent length 
(L), length of the femur (F), length of the tibia (T), length 
of the first toe (DP), and length of the metatarsal tubercle 
(CINT). These characteristics vary between the taxa and 
are traditionally used in the morphometry of water frogs 
(GUNTHER 1990, PLOTNER 2005). To minimize bias intro- 
duced by measurements, all frogs were measured by a sin- 
gle person (P. MIKULiCEK). The species status of every indi- 
vidual was determined by combining morphological char- 
acteristics and colour patterns in the field (SPEYBROECK 
et al. 2016). Additionally, the shape of the metatarsal tu- 
bercle was recorded, and individuals photographed. Only 
individuals with L > 41.50 mm (P. shqipericus), 45.20 mm 
(P. epeiroticus), 51.50 mm (P. ridibundus) and 52.24 mm 
(P. kurtmuelleri) were involved in the study. Males with this 
body length had clearly visible vocal sacs and were distin- 
guishable from females. 

A tissue sample (a drop of blood or toe clip) from each 
individual was fixed in 96% ethanol for molecular analy- 
sis. 


100 200 300 400km 


Genetic identification of Balkan species 


Species identification of Balkan individuals was based on 
the size polymorphism and sequence variation in serum al- 
bumin intron-1 (SAI-1) (PLOTNER et al. 2009, HAUSWALDT 
et al. 2012), and on microsatellite markers (see below). The 
size polymorphism of the species-specific SAI-1 fragments 
was examined on an agarose gel (HAUSWALDT et al. 2012). 
Identification of the Central-European P. ridibundus frogs 
was based on the combined molecular markers in our pre- 
vious works (MIKULÍČEK et al. 2014, 2015) and is not sub- 
jected to this study. 

Genomic DNA was extracted from blood or toe clips us- 
ing NucleoSpin® Tissue kit (Macherey-Nagel, Düren, Ger- 
many) following the manufacturer’s protocol. For ampli- 
fication of the SAI-1 fragment, the primers Pel-SA-Fi and 
Pel-SA-Ri were used (HAUSWALDT et al. 2012). PCR was 
carried out in a total volume of 10 ul with 5 ul of VWR Red 
Taq 2x mix (VWR, Radnor, PA, USA), 0.3 ul of each primer 
(10 uM), 3.4 ul of ddH,O, and 1 ul of DNA. PCR program 
was as follows: 2 min of initial denaturation at 94°C fol- 
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Figure 1. Localities of water frogs and species composition in the southwestern Balkans and Slovakia (A). The inset (B) shows sample 
sites from the southwestern Balkans in detail. 1 - Berat, 2 - Dushnik, 3 - Qazim Pali, 4 - Ksamil, 5 - Mbrostar, 6 - Divjakë - Ka- 
ravasté, 7 — Nishaj, 8 - Pocem, 9 - Xare, 10 - Zvërnec, 11 - Igoumenitsa, 12 - Ioannina, 13 - Virpazar, 14 - Brodské, 15 - Čunovo, 
16 - Devin, 17 - Sur, 18 - Šaštín-Stráže. 
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lowed by 35 cycles of denaturation for 30 s at 94°C, primer 
annealing at 57°C for 30 s, and elongation at 72°C for 1 min, 
with a final elongation at 72°C for 10 min (modified from 
Hauswa1pt et al. 2012). PCR products were run in agar- 
ose electrophoresis with a molecular-weight size marker. 
SAI-1 alleles of selected individuals were sequenced com- 
mercially in Macrogen Europe (Amsterdam, the Nether- 
lands), and their similarity with published sequences was 
compared using the program BLAST (https://blast.ncbi. 
nlm.nih.gov/). Obtained sequences are deposited in the 
NCBI GenBank under the accession numbers MW296101- 
105 for P. shqipericus; MW287331-335 for P. epeiroticus and 
finally, MW296106-111 for P. kurtmuelleri. Relationships 
among SAI-1 alleles were reconstructed using a parsimony 
network algorithm of TCS (CLEMENT et al. 2000) imple- 
mented in PopArt 1.7 (http:// http://popart.otago.ac.nz). 
All sequences were aligned and checked using Geneious 
Prime 2020 2.4 (Biomatters, Auckland, New Zealand). A 
final alignment involved 45 sequences of the size of 597 
bp. Pelophylax perezi (LOPEZ-SEOANE, 1885) was used as 
an outgroup based on previous studies (LYMBERAKIS et al. 
2007). 

Microsatellites were amplified in multiplex PCRs. Most 
of individuals were analyzed for five microsatellite loci: 
Ga1a19 (CHRISTIANSEN 2009), RICA1b5 (GARNER et al. 
2000), Rrido13A (Horz et al. 2001), RICA2a34 and Rrid- 
135A (CHRISTIANSEN & REYER 2009). PCR reactions were 
performed in a total volume of 10 ul and consist of 5 ul of 
Qiagen Microsatellite PCR Master mix (Qiagen, Hilden, 
Germany), 0.1 ul of each primer (10 uM), 3 ul of ddH,O 
and 1 ul of DNA. In the case of 31 individuals of P. shqiperi- 
cus and 40 individuals of P. kurtmuelleri, 13 microsatellite 
loci were used for species identification: Ga1ai9 (CHRIS- 
TIANSEN 2009), RICA1bs, RICA2a34F, RICA5 (GARNER et 
al. 2000), RICA1b6, ReiCagaio, Re2Caga3 (ARIOLI et al. 
2010), Resi4, Res22 (ZEISSET et al. 2000), Rrido13A, Rrid- 
082A (Horz et al. 2001) and Rridi69A. PCR program was 
modified from CHRISTIANSEN & REYER (2009): 5 min of 
initial denaturation at 95°C followed by 30 cycles of dena- 
turation for 30 s at 95°C, 60°C for 90 s and 72°C for 1 min, 
with a final extension at 60°C for 30 min. Microsatellite 
fragments were run on an automated ABI genetic analyser. 

The Bayesian clustering implemented in the program 
Structure 2.3.4 (PRITCHARD et al. 2000) was used for mic- 
rosatellite data analyses. In the first analysis, all three spe- 
cies were analysed together for five microsatellite loci. Sub- 
sequently, the frogs were analysed in two datasets involv- 
ing sympatrically living species pairs P. shgipericus/P. kurt- 
muelleri (the first dataset based on thirteen loci), and 
P. epeiroticus/P. kurtmuelleri (the second dataset based on 
five loci), respectively. The admixture and non-correlated 
allele model was used to calculate the parameter Q, that 
is, the proportion of an individual’s genome originating in 
one of the inferred clusters (K), corresponding to the spe- 
cies P. shqipericus, P. epeiroticus and P. kurtmuelleri. Fixed 
K = 3 was chosen when all three species were analysed to- 
gether, fixed K = 2 when species pairs were analyzed. The 
analyses were based on runs of 10° iterations, following a 


108 


burn-in period of 10° iterations. Five independent runs 
were carried out to test the accuracy of the results. Besides 
Bayesian analyses, the Discriminant Analysis of Principal 
Components (DAPC) based on five microsatellite loci was 
performed. The ‘adegenet' package (JOMBART et al. 2020) 
implemented in the R statistical environment (R Core 
Team 2019) was used for this analysis (for more details 
about the DAPC see below). 


Morphometric analyses 


Statistical analysis of morphometric characteristics was 
performed in the R statistical environment (R Core Team 
2019). A suitable data format was created using 'tidyverse' 
(WICKHAM & RSTUDIO 2019) and ‘janitor’ (FIRKE et al. 
2019) packages. Seven morphometric indices (L/F, L/T, 
L/DP, L/CINT, F/T, T/CINT, DP/CINT) were calculated 
(GUNTHER 1990, PLOTNER 2005). Before conducting anal- 
ysis, correlation, outliers and the assumption of normal- 
ity were checked via the 'PerformanceAnalytics' package 
(PETERSON et al. 2020). Two-way ANOVA was performed 
to examine differences in calculated proportions among 
species and sexes. When ANOVA was significant, Tukey 
HSD post-hoc tests were used to perform multiple pair- 
wise comparisons among the species. Both analyses were 
performed with respect to 'rstatix’ package syntax (Kas- 
SAMBARA 2020). Identification of clusters and distances in 
morphospace was performed via the Discriminant Anal- 
ysis of Principal Components (DAPC) from ‘adegenet' 
package (JOMBART et al. 2020). This method finds the lin- 
ear combinations of morphometric variables that have the 
largest between-group variance and the smallest within- 
group variance. DAPC uses the transformation of varia- 
bles by the Principal Component Analysis (PCA) prior to 
performing the Discriminant Analysis (DA). This ensures 
that variables are uncorrelated, thus making it a suitable 
methodological approach for morphological data that usu- 
ally demonstrate a high level of correlation. Based on the 
results of DAPC, group membership probabilities for each 
individual were computed. All visualizations were created 
using the 'ggplot2' (WicKHaM et al. 2019), 'ggpubr' (Kas- 
SAMBARA 2019), 'adegenet' (JOMBART et al. 2020) and 'grid- 
Extra’ (AUGUIE & ANTONOV 2017) packages. 


Results 
Genetic identification of Balkan species 


Amplification of the SAI-1 revealed three PCR fragments 
of the length ~720 bp, ~840 bp, and ~1 200 bp, which cor- 
responded to the species P. kurtmuelleri, P. epeiroticus, and 
P. shqipericus, respectively. The BLAST supported a correct 
species-specificity of sequenced PCR fragments. SAI-1 al- 
leles formed species-specific clusters in a parsimony net- 
work (Supplementary Figure S1). Four individuals were 
heterozygous in SAI-1 (Supplementary Table S1). A female 
captured near Ksamil, Albania (ID: MIK3233), and a male 
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from the Ioannina Lake, Greece (MIK3309) carried one al- 
lele specific for P. kurtmuelleri and one for P. epeiroticus. 
A male from Divjakë - Karavasté, Albania (MIK3160) had 
one P. kurtmuelleri allele and one P. shqipericus allele. Fi- 
nally, a female captured in Pocem, Albania (JAB5636), 
outside from P epeiroticus distribution range, carried 
two fragments indicating it was a hybrid between P kurt- 
muelleri and P. epeiroticus. In this case, the sequence of the 
fragment ~840 bp (MW296111) was not homologous with 
P. epeiroticus but matched P. ridibundus individuals from 
Kamchatka, Russia (99.50% identity, accession number 
KX503320.1), Aliartos, Greece (99.37% identity, accession 
number MF667644.1), and Atyrau, Kazakhstan (99.37% 
identity, accession number HE858213.1). This SAI-1 allele 
was assigned to the P ridibundus cluster in a parsimony 
network (Supplementary Figure Si), and thus a female 
JABs5636 was not considered as a hybrid but as a P. kurt- 
muelleri individual with a ridibundus-specific allele 

The program Structure assigned individuals to the clus- 
ters corresponding to the species P. shqipericus, P. epeiro- 
ticus and P. kurtmuelleri (Supplementary Table S2). The 
following four individuals revealed an admixed genome 
and were excluded from morphological analyses. A male 
from the Ioannina Lake (MIK3309) carried 57.5% of ge- 
nome from P. kurtmuelleri and 42.5% from P. epeiroticus, 
and was considered an Fi hybrid. Two individuals from 
Ksamil (MIK3233) and Igoumenitsa (MIK3277) carried 
75.0% and 77.5% of genome from P. kurtmuelleri, the re- 
maining allelic variation was derived from P. epeiroticus, 
and both were considered backcrossed hybrids. However, 
MIK3277 was assigned to the P. kurtmuelleri cluster when 
all three species were analyzed together (with fixed K = 3; 
Supplementary Table S2). A male from Divjakë - Karavasté 
(MIK3160) possessed 30.1% of genome from P. kurtmuelleri 
and 69.9% from P. shqipericus, and was also evaluated as 
a backcrossed hybrid. A female from Poçem (JAB5636), 
which carried a SAI-1 allele specific for P. ridibundus, was 
clearly assigned to the P. kurtmuelleri cluster. 

Microsatellite-based DAPC revealed three spatially sep- 
arated clusters corresponding to particular species. All hy- 
brid individuals but MIK3277 showed intermediate posi- 
tion in multivariate space (Fig. 3B). 


Morphometric analyses 


Morphometric characteristics and indices for both males 
and females of each species are summarized in Tables 2 and 
3, and in Figure 2. The ANOVA revealed significant differ- 
ences between P epeiroticus and P. kurtmuelleri in six out 
of seven morphometric indices in males and in all seven in 
females (Table 4). Pelophylax shqipericus and P. kurtmuel- 
leri differed significantly from each other at six indices for 
both sexes. Pelophylax kurtmuelleri and P. ridibundus dif- 
fered significantly at six indices in males, whereas only four 
indices were significantly different in females. Interesting- 
ly, P. epeiroticus and P. shqipericus were the most similar 
species pair. No significantly different indices were found 


in males, and only two significantly different indices were 
in females. The significantly different indices across species 
were L/T, L/CINT, and DP/CINT (Tab. 4). 

Four individuals identified as nuclear hybrids shared 
morphometric indices with both parental taxa (Supple- 
mentary Table S1). A specimen MIK3309 possessing 57.5% 
of genome from P. kurtmuelleri and 42.5% from P. epeiroti- 
cus had three and four indices similar to their parents, re- 
spectively. Specimens MIK3233 and MIK3277 (with about 
three-quarters of genome from P kurtmuelleri) had the 
most indices similar to P. kurtmuelleri. Finally, a hybrid, 
MIK3160, composed of 69.9% from P. shqipericus genome 
and of 30.1% from P. kurtmuelleri had three and four indi- 
ces similar to their parents, respectively. 

Three principal components were retained for DAPC 
analysis. The first PC had an eigenvalue of 12.7, accounting 
for 95% of the variance (Supplementary Table S3). The oth- 
er two retained components (PC2, PC3) had significantly 
lower eigenvalues (Supplementary Table S3). Overall, re- 
tained variance by all three principal components was 99%. 
The first principal component had noticeable loadings on 
L/CINT and T/CINT variables. The second principal com- 
ponent displayed strong loading on L/DP variable, and the 
third retained PC had strong loading on T/CINT variable 
(Supplementary Table S3). The examination of morphos- 
pace created via PCA (Supplementary Figure S2) showed 
a high overlap of all studied species. Similarly, the applica- 
tion of Discriminant Analysis (DA) introduced a consider- 
able amount of overlap in all species (Fig. 3A). Both cen- 
troids of P epeiroticus and P. shqipericus and P. ridibundus, 
P. kurtmuelleri appear to be closely situated in morpho- 
space. Both P. epeiroticus and P. shqipericus, however, ap- 
pear to be more distant in morphospace from P. kurt- 
muelleri. The hybrids either fell to the centroids of the pa- 
rental species (MIK3160, MIK3277) or showed an interme- 
diate position in multivariate morphospace (MIK3233 and 
MIK3309). 


Colouration and qualitative morphological traits 


The dorsal colouration of P. shqipericus (Fig. 4) varied from 
bright green to olive-brown or even dark brown. The dor- 
sal pattern consisted of relatively large irregular dark green, 
olive, or dark grey spots. Almost uniformly coloured indi- 
viduals with a strongly reduced spot pattern were rare (ap- 
prox. 3% of individuals). A mid-dorsal yellow or light green 
stripe was present in the majority of individuals. Tympana 
can be green, olive, bronze, or dark brown with a lighter 
centre and darker surroundings or unicolour. Dorsolateral 
folds were visible but low, often distinguishable by a dis- 
tinct colour. Male vocal sacs are light grey or olive green 
(Fig. 7A). Lateral sides are usually paler, marked with dark- 
er spots with diffused edges and yellow pigment. The up- 
per side of the legs is coloured and spotted similar to the 
dorsum. Several darker transversal spots were present in 
the anterior part of the thighs, dark grey to black merged 
spots with yellow pigment were typical for a posterior part 
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Table 2. Summary statistics of morphometric characters measured in three water frog species (genus Pelophylax) from the southwestern 
Balkans and P ridibundus from Central Europe. Min - minimal value, Max - maximal value, SD - standard deviation, SE - standard error. 


Character P. epeiroticus P. kurtmuelleri P ridibundus P. shqipericus 
Males Females Males Females Males Females Males Females 
N=25 N=26 N=57 N=43 N=57 N=32 N=36 N=24 
L 
Mean 55.70 70.06 67.69 71.67 80.48 85.20 52.63 54.27 
Min 45.20 46.87 52.24 54.69 51.50 50.00 41.50 42.00 
Max 72.10 82.50 85.10 104.20 111.40 117.30 66.51 81.00 
Median 53.90 71.35 67.50 71.80 86.40 89.40 51.45 52.45 
SD 7:12 8.36 9.42 10.81 15.06 21.76 7.03 8.87 
SE 1.42 1.64 1.25 1.65 1.99 3.85 1:17 1.81 
F 
Mean 27.22 33.10 34.23 36.12 40.21 42.16 25.66 25.63 
Min 22.00 22.15 26.59 28.33 24.95 22.60 20.00 19.90 
Max 37.27 40.20 44.30 50.50 53.10 59.65 32.70 39.30 
Median 26.28 33.05 33.86 35.42 43.45 44.05 25.25 25.00 
SD 3.70 4.16 4.78 4.93 7.50 11.16 3.55 4.19 
SE 0.74 0.82 0.63 0.75 0.99 1.97 0.59 0.86 
T 
Mean 28.49 34.23 37.51 39.95 42.90 45.23 27.04 27.07 
Min 23.20 22.93 28.69 30.95 26.70 21.90 20.60 21.60 
Max 36.31 39.20 48.00 56.00 56.35 62.50 34.30 41.00 
Median 27.60 35.18 37.11 39.36 46.85 48.91 26.85 26.00 
SD 3.25 3.92 4.91 5.46 8.19 12.50 3.75 4.33 
SE 0.65 0:77 0.65 0.83 1.08 2.21 0.63 0.88 
DP 
Mean 8.86 10.81 9.51 10.06 11.49 12.20 8.01 8.14 
Min 6.44 7.31 6.92 Asbo 6.20 5.30 6.20 6.00 
Max 11.72 12.70 13.10 14.60 15.35 16.70 10.40 11.50 
Median 8.75 10.88 9.46 10.00 12.25 13.80 8.20 8.05 
SD 1.35 1.46 1.39 1.45 2.36 3.53 1.05 1.24 
SE 0.27 0.29 0.18 0.22 0.31 0.62 0.18 0.25 
CINT 
Mean 2.84 2.90 3.67 3.86 3.95 4.14 2.57 2.49 
Min 2.10 1.90 2.40 2.74 2.30 2.05 1.80 1.90 
Max 3.98 3.97 5.59 5.28 3:99 6.87 4.12 3.70 
Median 2:72 2.90 3.50 3.81 3.85 4.35 2.50 2.40 
SD 0.49 0.48 0.61 0.67 0.81 1.00 0.59 0.51 
SE 0.10 0.09 0.08 0.10 0.11 0.18 0.10 0.10 


(Fig. 8A). The venter was usually white or yellowish-white 
with diffused darker spots in some individuals and with 
a higher amount of yellow pigment. Similar colouration 
was also present on the ventral side of the thighs (Fig 9A). 
Rarely, there was brown spots in the anterior part of the 
venter. The metatarsal tubercle was either triangular or 
oblique, usually low, but sometimes higher and reached a 
maximum height in its distal part. The yellow or orange 
pigment was present in the foot webbing (Fig. 7A). 

Dorsal colouration of P. epeiroticus (Fig. 5) formed a 
scale from yellowish green to dark green or olive green 
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with relatively large, irregular, dark green to black spots. 
However, individuals without spots were also recorded. 
A yellow or light green mid-dorsal stripe was present in 
the majority of individuals. Dorsolateral glandular folds 
were low, but clearly visible and often olive or brown. Al- 
though, these folds can bear the same colour as the sur- 
rounding skin. The tympanum was either green, bronze, 
or brown, often with a darker margin and lighter cen- 
tre. Vocal sacs of males were dark grey to almost black 
(Fig. 7B). Lateral sides were usually paler than, or simi- 
larly coloured to, a dorsal side with the yellow pigment. 


Morphological differentiation of Pelophylax from the southwestern Balkans 


Table 3. Summary statistics of morphometric indices in three water frog species (genus Pelophylax) from the southwestern Balkans 
and P. ridibundus from Central Europe. Min - minimal value, Max - maximal value, SD - standard deviation, SE - standard error. 


Index P. epeiroticus P. kurtmuelleri P. ridibundus P. shqipericus 
Males Females Males Females Males Females Males Females 
N=25 N=26 N=57 N=43 N=57 N=32 N=36 N=24 


L/F 
Mean 2.05 2.12 1.98 1.98 2.00 2.03 2.05 212: 
Min 1.91 1.99 1.80 1.82 1.88 1.86 1.89 1.94 
Max 2.19 2.29 2.17 2.00 2.21 2.21 2.21 2.28 
Median 2.05 2.11 1.97 2.00 2.00 2.02 2.04 2.12 
SD 0.07 0.08 0.07 0.07 0.08 0.08 0.07 0.08 
SE 0.01 0.02 0.01 0.01 0.01 0.01 0.01 0.02 
L/T 
Mean 1.95 2.05 1.80 1.79 1.88 1.90 1.95 2.00 
Min 1.86 1.91 1.68 1.67 1.71 1.75 1.85 1.84 
Max 2.03 2.12 1.93 1.96 2.10 2.28 2.14 2.11 
Median 1.96 2.05 1.81 1.80 1.87 1.88 1.94 1.99 
SD 0.05 0.05 0.06 0.06 0.09 0.12 0.06 0.07 
SE 0.01 0.01 0.01 0.01 0.01 0.02 0.01 0.01 
L/DP 
Mean 6.32 6.50 7.14 7.13 7.07 7.15 6. 6.67 
Min 5.57 6.04 6.01 5.78 5.89 5.99 5.77 5.89 
Max 7.39 7.05 7.91 8.10 8.61 9.47 7.27 7.40 
Median 6.33 6.43 7.16 7.20 7.02 7.10 6.52 6.64 
SD 0.40 0.27 0.46 0.49 0.62 0.89 0.40 0:37 
SE 0.08 0.05 0.06 0.07 0.08 0.16 0.07 0.08 
L/CINT 
Mean 19.83 24.53 18.64 18.74 20.64 20.77 21.00 22.08 
Min 15.82 19.32 14.76 14.98 15.14 14.14 15.36 17:71 
Max 24.62 32.74 24.83 23.67 28.66 28.15 28.17 25.45 
Median 19.64 24.42 18.28 18.71 20.21 21.48 21.31 22.34 
SD 2.09 3.20 2.26 2.15 3.07 3.59 2.90 2.36 
SE 0.42 0.63 0.30 0.33 0.41 0.64 0.48 0.48 
F/T 
Mean 0.95 0.97 0.91 0.90 0.94 0.94 0.95 0.95 
Min 0.89 0.92 0.85 0.86 0.86 0.86 0.89 0.90 
Max 1.03 1.04 0.97 0.96 1.03 1.04 1.00 1.00 
Median 0.95 0.96 0.91 0.90 0.94 0.93 0.95 0.94 
SD 0.04 0.03 0.03 0.02 0.04 0.04 0.02 0.02 
SE 0.01 0.01 0.00 0.00 0.01 0.01 0.00 0.00 
T/CINT 
Mean 10.16 11.99 10.34 10.46 11.03 11.02 10.77 11.02 
Min 8.12 9.35 7.82 8.15 7.63 6.81 8.19 8.78 
Max 13.02 15.65 13.40 12.48 15.46 15.51 14.44 13.33 
Median 10.12 11.97 10.33 10.36 10.94 11.20 10.78 11.14 
SD 1.17 1.50 1.22 1.16 1.92 2.32 1.41 1.22 
SE 0.23 0.29 0.16 0.18 0.25 0.41 0.24 0.25 
DP/CINT 
Mean 3:15 3.78 2.63 2.64 2.96 2.98 Dal 3.32 
Min 2.55 2.87 1.97 2.04 1.77 1.65 2.19 2.44 
Max 4.03 5.13 3.53 3.48 4.62 4.43 4.56 4.24 
Median 3.11 3.63 2.60 2.57 2.86 3.03 3.20 3.30 
SD 0.39 0.49 0.38 0.39 0.62 0.76 0.50 0.41 
SE 0.08 0.10 0.05 0.06 0.08 0.13 0.08 0.08 
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Figure 2. Comparison of morphometric indices among three Balkan water frog species (genus Pelophylax) and central-European 
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P. ridibundus. * - p < 0.05, ** - p < 0.005, *** - p < 0.0005, ns — not significant. 
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Morphological differentiation of Pelophylax from the southwestern Balkans 


Table 4. Differences between water frog species (genus Pelophylax) in morphometric indices revealed by ANOVA and Tukey HSD post- 
hoc tests with statistically significant values in bold. S - Significance, * - p < 0.05, ** - p < 0.005, *** - p < 0.0005, ns — not significant. 


Species 1 Species 2 Sex Variable p value S Species 1 Species 2 Sex Variable p value S 

P. epeiroticus P shqipericus Males L/F 1.00 e* ns P. epeiroticus P. shqipericus Females L/F 1.00 e*° ns 
P. epeiroticus P. shqipericus Males L/T 2.99 e! ns P. epeiroticus P. shqipericus Females L/T 2.71 e1 ns 
P. epeiroticus P shqipericus Males  L/DP 1.61 e? ns P. epeiroticus P. shqipericus Females L/DP 7.00 e? ns 
P. epeiroticus P. shqipericus Males L/CINT 5.37e7 ns P. epeiroticus P. shqipericus Females L/CINT 1.48 e° * 

P. epeiroticus P shqipericus Males F/T 2.03 e! ns P. epeiroticus P. shqipericus Females F/T 1.26 e? ns 
P. epeiroticus P. shqipericus Males T/CINT 9.34e' ns P. epeiroticus P. shqipericus Females T/CINT 1.55 e? ns 
P. epeiroticus P. shqipericus Males DP/CINT 1.79 e? ns P. epeiroticus P. shqipericus Females DP/CINT 1.49 e? * 

P. epeiroticus P kurtmuelleri Males L/F 1.36 e? *** P. epeiroticus P. kurtmuelleri Females L/F 529er 0A 
P. epeiroticus P kurtmuelleri Males L/T 6.68 e1 *** P. epeiroticus P. kurtmuelleri Females L/T 8.66 e17 *** 
P. epeiroticus P kurtmuelleri Males L/DP 9.09 e1 *** P. epeiroticus P. kurtmuelleri Females L/DP 1.05 e4 *** 
P. epeiroticus P. kurtmuelleri Males L/CINT 9.26 e™ *** P. epeiroticus P. kurtmuelleri Females L/CINT 1.91 e”? *** 
P. epeiroticus P kurtmuelleri Males F/T 7.04 e” *** P. epeiroticus P. kurtmuelleri Females F/T 7.50 e1 *** 
P. epeiroticus P. kurtmuelleri Males T/CINT 6.92 e? ns P. epeiroticus P. kurtmuelleri Females T/CINT 1.25 e° ** 
P. epeiroticus P. kurtmuelleri Males DP/CINT 7.06 e” *** P. epeiroticus P. kurtmuelleri Females DP/CINT 2.70 e™ *** 
P. kurtmuelleri P. shqipericus Males L/F 7.89 e1 +* P. kurtmuelleri P. shqipericus Females L/F 1.64 e°? *** 
P. kurtmuelleri P. shqipericus Males L/T 6.68 e1 *** P. kurtmuelleri P. shqipericus Females L/T 4.15 e74 *** 
P. kurtmuelleri P. shqipericus Males L/DP 3.88 e8 *** P. kurtmuelleri P. shqipericus Females L/DP 1.09 e? * 

P. kurtmuelleri P. shqipericus Males L/CINT 4.26e8 *** P. kurtmuelleri P. shqipericus Females L/CINT 6.23 e5 *** 
P. kurtmuelleri P. shqipericus Males F/T 4.59 e1 *** P. kurtmuelleri P. shqipericus Females F/T 3.57 e6 *** 
P. kurtmuelleri P. shqipericus Males T/CINT 2.33e7 ns P. kurtmuelleri P. shqipericus Females T/CINT 5.24e? ns 
P. kurtmuelleri P. shqipericus Males DP/CINT 1.33 e™ *** P. kurtmuelleri P. shqipericus Females DP/CINT 1.25 e° *** 
P. kurtmuelleri P. ridibundus Males L/F 9.15 e° ** P. kurtmuelleri P. ridibundus Females L/F 5.39 e? ns 
P. kurtmuelleri P. ridibundus Males L/T yae o P. kurtmuelleri P. ridibundus Females L/T 2.07 e7 *** 
P. kurtmuelleri P. ridibundus Males L/DP 9.59 e! ns P. kurtmuelleri P. ridibundus Females L/DP 1.00 e*° ns 
P. kurtmuelleri P. ridibundus Males L/CINT 5.71e6 *** P. kurtmuelleri P. ridibundus Females L/CINT 1.49 e? * 

P. kurtmuelleri P. ridibundus Males F/T 5.87 e° *** P. kurtmuelleri P. ridibundus Females F/T 8.32 e5 *** 
P. kurtmuelleri P. ridibundus Males T/CINT 2.18e? * P. kurtmuelleri P. ridibundus Females T/CINT 4.51 e7 ns 
P. kurtmuelleri P. ridibundus Males DP/CINT 3.73 e5 *** P. kurtmuelleri P. ridibundus Females DP/CINT 3.50 e? * 


Lateral spots were usually present, dark grey or black, 
always with diffused edges. Yellow pigment is present 
on the margins of the ventral side. Legs were dorsally 
green or yellowish with olive, dark grey, or black irregu- 
lar spots. These spots were transversal and merged or are 
separated (Fig. 8B). The ventral side was uniformly white 
to creamy, sometimes with grey or brown marbling. The 
ventral part of the thighs was uniformly white to creamy 
or covered by marbling (Fig. 9B). The metatarsal tuber- 
cle was predominantly very small and low and either the 
shape of a scalene triangle or oblique. Orange or yellow 
pigment was characteristically present in the foot web- 
bing (Fig. 7B). 

In P. kurtmuelleri (Fig. 6), the dorsal part of the body 
was green to olive green or brown. Green colouration was 
usually situated in the anterior part of the body; the poste- 
rior part was usually brown or grey. Brown, grey or black, 
mostly irregular spots were distributed over the dorsal 
part, uniformly coloured individuals were rare. A yellow 
or light green mid-dorsal stripe was present in the major- 


ity of individuals. Dorsolateral folds were clearly visible, 
often distinguishable from the surrounding skin by a dis- 
tinct colour. Tympana were olive, bronze, brown, or rare- 
ly green, with a paler centre. Vocal sacs in males were ol- 
ive, grey, or dark grey (Fig. 7C). Lateral sides have a lighter 
tint compared to the dorsum and usually many irregular 
solitaries or merged black spots with sharp (not diffused) 
edges. Hind legs were usually olive or brown with dark- 
er transversal spots in the anterior part and merged dots 
in the posterior part, forming a web-like pattern. Yellow 
pigment was absent (Fig. 8C). The ventral side was mainly 
white without any pattern, rarely with brown or dark grey 
spots situated between the forelegs, in the mandibular area, 
or across the whole belly. The ventral side of the legs was 
predominantly white and spotless (Fig. 9C). The metatarsal 
tubercle was low but higher than in P. epeiroticus, and ei- 
ther triangular, or oblique. The foot webbing lacks any yel- 
low or orange pigmentation (Fig. 7C). Table 5 characterizes 
the Balkans species living in sympatry for their identifica- 
tion in the field. 
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Discussion 


Proper taxon identification is a prerequisite of any biologi- 
cal study, with morphological and morphometric charac- 
ters being the first to be examined. Univariate (ANOVA) 
and multivariate (DAPC) analyses of water frogs from 
the southwestern Balkans revealed substantial differenc- 
es between the sympatric species pairs, P kurtmuelleri/ 
P. shqipericus, and P. kurtmuelleri/P. epeiroticus. The most 
useful characteristics for the species identification are mor- 
phometric indices L/CINT, T/CINT, L/T, and DP/CINT. 
Key morphological characters for identifying P. epeiroticus 
are the presence of yellow pigment in the flanks, venter, 
groins and thighs; lateral spots absent or with diffused edg- 
es and a small metatarsal tubercle. Similar characters also 
fit to P shqipericus, with exception of the metatarsal tuber- 
cle, which is bigger and higher compared to P. kurtmuel- 
leri. On the other hand, yellow pigment in P. kurtmuelleri 
is usually absent; lateral spots with sharp edges are present 
and finally, the metatarsal tubercle is bigger compared to 
P. epeiroticus, but smaller than in P. shqipericus. Such data 
extends the previous allozyme and immunological mark- 
ers (HOTZ & UZZELL 1982, Hotz et al. 1987), or bioacous- 
tic parameters (species-specific differences in mating calls; 
SCHNEIDER et al. 1984, 1993) for species identification. The 
advantage of our morphological and morphometric indi- 


DA eigenvalues 


+ Pelophylax epeiroticus 
* Pelophylax kurtmuelleri 
+ Pelophylax ridibundus 

+ Pelophylax shqipericus 


ces is their applicability where no molecular analyses are 
available, allowing for non-invasive field identification. 
The availability of identification markers may significantly 
help in potential sympatric populations to revise species 
distribution ranges and aid in species conservation. 


Validity of morphometric and morphological tools 
to differentiate water frog species 


Our study on a number of species across a wide distribution 
range confirmed a general validity of morphometric indices 
obtained by Hotz & UZzELL (1982) and Horz et al. (1987) 
from Montenegro and SCHNEIDER et al. (1984, 1993) from 
Greece. The only exceptions were indices, where CINT was 
involved. Our more extensive sampling shows that the val- 
ues T/CINT and L/CINT in P. shqipericus and P. kurtmuel- 
leri, and T/CINT and DP/CINT in P epeiroticus were high- 
er than the values presented in the aforementioned studies. 

According to ANOVA, the most significantly different 
indices between the species pairs were L/T, L/CINT, and 
DP/CINT for P. epeiroticus and P. kurtmuelleri, as also pro- 
posed by SCHNEIDER et al. (1984), and L/T, L/F, F/T and 
DP/CINT for P. shqipericus and P. kurtmuelleri. ANOVA 
results are corroborated by multivariate analyses. In gen- 
eral, a multivariate approach revealed considerable overlap 
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Figure 3. Discriminant analysis of principal components (DAPC) based on morphological (A) and microsatellite markers (B). Both 
analyses were carried out for all three (microsatellites) or four (morphological markers) species (upper panels) and subsequently for 
sympatrically living species pairs (lower panels). The position of hybrid individuals is displayed. A specimen MIK_3277 was identified 


as a hybrid only in a Bayesian analysis (see Results). 
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Figure 4. Variability of colour pattern in Pelophylax shqipericus from Nishaj, Albania (A, B, C); Velipojë, Albania (D, E, Ẹ); Virpazar, 
Montenegro (G) and Divjaké, Albania (H). 
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Figure 5. Variability of colour pattern in Pelophylax epeiroticus from Igoumenitsa, Greece (A, B); Ioannina, Greece (C, D, E); Corfu, 
Greece (F); Kalogria, Greece (G) and Stjar, Albania (H). 
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Figure 6. Variability of colour pattern in Pelophylax kurtmuelleri from Pogem, Albania (A, B); Qazim Pali, Albania (C, D); Shkalla, 
Albania (E); Mbrostar, Albania (F), Velipojë, Albania (G) and Perbreg, Albania (H). 
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among all species after the analysis of principal components 
(PCA). However, after the following DAPC, syntopically 
living P. kurtmuelleri/P. shqipericus or P. epeiroticus seem to 
be more separated in morphospace than other species pairs. 
The morphological differentiation between the Europe- 
an water frog species is relatively narrow despite their deep 
phylogenetic divergence and high genetic differentiation 
as we also show from the studied regions. The valued dis- 
crimination characters are differences in the length of the 
body, hind legs, the first toe, and the metatarsal tubercle, as 
well as the head shape, the shape of the metatarsal tubercle 
and colouration (GUNTHER 1990, PLOTNER 2005). Thus, 
the water frog species are rather challenging to correctly 
identify in the field, and this holds for the southwestern 
Balkans as well (e.g. SzABOLCs et al. 2017). So, how to dif- 
ferentiate between sympatric P. kurtmuelleri/P. epeiroticus, 
and P. kurtmuelleri/P. shqipericus individuals correctly? 


Species and hybrid identification in the field 


According to our observations, P. epeiroticus differs from 
P. kurtmuelleri (Supplementary Figure S3) in the follow- 
ing traits: the presence of the yellow pigment in the groin, 
flanks, belly edges and upper as well as lower parts of the 


thighs, lateral spots are absent or, if present, with diffused 
edges, yellow or orange feet webbing, and a very small 
and oblique metatarsal tubercle. Moreover, vocal sacs of 
P. epeiroticus males are usually darker when compared to 
P. kurtmuelleri. 

In morphometry, the best indices separating these two 
species are L/T (P. epeiroticus > 2, P. kurtmuelleri < 1.9), 
L/CINT (P epeiroticus > 23.7, P. kurtmuelleri < 15.8), and 
DP/CINT (P epeiroticus > 3.5, P. kurtmuelleri < 2.9). These 
findings agree with other authors (HOTZ & UzzELL 1982, 
SCHNEIDER et al. 1984, SPEYBROECK et al. 2016). Addition- 
ally, P epeiroticus should have a shorter tibia and femur 
(Hotz & UZZELL 1982), which is also evident in our data. 
On the other hand, we cannot confirm the unique green 
tympana colouration suggested by SCHNEIDER et al. (1984) 
as our observations also included P. epeiroticus individuals 
with bronze tympana, or even a combination of both. 

Comparing another syntopically living species pair, 
P. shqipericus and P. kurtmuelleri, the best discriminative 
characteristics for P. shqipericus are yellow pigmentation 
in the groin, flanks, belly edges and on the thighs, lateral 
spots with diffused edges or absent, pronounced yellow or 
orange colouration of feet webbing, a relatively high and 
usually triangular metatarsal tubercle, and light greyish or 
greenish colouration of vocal sacs of males. The best indi- 


Figure 7. The shape of the metatarsal tubercle, vocal sacs colouration and the colouration of the foot webbing of Pelophylax shqipericus 
(A), P. epeiroticus (B) and P. kurtmuelleri (C). 
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Figure 8. A pattern of lateral spots and yellow pigment in venter margins and upper thighs of Pelophylax shqipericus (A), P. epeiroticus 
(B) and P. kurtmuelleri (C). 


Figure 9. A ventral pattern and colouration of Pelophylax shqipericus (A), P. epeiroticus (B) and P. kurtmuelleri (C). 


119 


PETR PAPEZiK et al. 


Table 5. Comparison of morphometric indices and qualitative morphological traits between sympatrically living species of water-frogs 
(genus Pelophylax) from the southwestern Balkans. ' - external morphological characteristics according to SPEYBROECK et al. (2016). 


Characteristics P. epeiroticus P. kurtmuelleri P. shqipericus P. kurtmuelleri 

L/T > 2.0 < 1.9 > 2.0 < 18 

L/CINT > 23:7 < 15.8 > 23.7 < 17.7 

T/CINT overlapped values overlapped values > 12.5 < 88 

DP/CINT > 3.5 < 2.9 > 3.5 < 2.4 

Ventral colouration Uniformly white or Uniformly white or with Usually white or Uniformly white or with 
(Fig. 8) creamy or with dark dark spots, which can be yellowish-white, rarely dark spots, which can be 


marbling 


or in the lower jaw and 


thighs 
Yellow pigment 


(Fig. 7) venter, upper and lower upper thighs 
thighs 

Lateral spots Absent or with diffused 

(Fig. 7) edges 


Shape and size of the 
metatarsal tubercle 


The shape of scalene 
triangle or oblique; 


situated over the belly 


with diffused dark spots 


situated over the belly 


Present in flanks, groins, Usually absent, rarely in 


Present, with sharp edges Absent or with diffused 


Triangular or oblique; 
low and small, but higher oblique; higher 


or in the lower jaw and 
thighs 
Present in flanks, groins, Usually absent, rarely in 
venter, upper and lower upper thighs 
thighs 
Present, with sharp edges 
edges 
Mostly triangular, rarely Triangular or oblique; 
lower 


Yellow to orange Olive, brown or dark 


(Fig. 6) usually very low and and larger than in 
small P. epeiroticus 

Foot webbing Yellow to orange Olive, brown or dark 

colouration (Fig. 6) grey 


Vocal sacs colour Dark grey, dark olive, 


(Fig. 6) black 

Darker than in 

P. kurtmuelleri 
Head shape’ Medium length and blunt Long with a sharp snout 
Skin! Smooth Rough with warts and 


ridges 
Heel joint when leg 
pulled forwards! 


Reaches the eyes 


Thigh colour during 
breeding period’ 


Yellow-spotted 
rarely yellow 


Olive, grey, dark grey 


Grey or white-spotted, 


grey 
Light grey or olive green Olive, grey, dark grey 
Lighter than in 
P. kurtmuelleri 


Medium length and sharp Long with sharp snout 
Smooth Rough with warts and 


ridges 


Reaches beyond the eyes, Reaches beyond the eyes, Reaches beyond the eyes, 
usually beyond the snout usually beyond the snout usually beyond the snout 


Yellow-spotted Grey or white-spotted, 


rarely yellow 


ces to distinguish between the species are L/T (P. shqiperi- 
cus > 2.0, P. kurtmuelleri < 1.8), L/CINT (P. shqipericus 
> 23.7, P. kurtmuelleri < 17.7), T/CINT (P shqipericus > 12.5, 
P. kurtmuelleri < 8.8) and DP/CINT (P. epeiroticus > 3.5, 
P kurtmuelleri < 2.4). Diffused lateral spots in P. shqipericus 
can continue to the ventral side and form dark marbling in 
contrast to P. kurtmuelleri, which has dark solid spots usu- 
ally only on the jaw or thigh area, whereas the venter is uni- 
formly white. SCHNEIDER et al. (1984), Horz et al. (1987) 
and SPEYBROECK et al. (2016) reported other differences 
between the Balkan water frog species like an extension of 
the webbing to toe tips with an indentation between, the 
smoothness of the skin, shape of the head or the position 
of the heel joint when the leg is pulled forward, but none of 
these were measured in our study. 

The morphological differentiation between P. ridibun- 
dus and P. kurtmuelleri were less pronounced. All seven 
measured morphometric indices to differentiate P kurt- 
muelleri from the Central European P. ridibundus are cryp- 
tic. However, we cannot exclude that a comprehensive 
morphological analysis including more external, skeletal 
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and other phenotypic characteristics would reveal a cor- 
relation between phylogenetic and morphological diver- 
gence and would shed light upon the taxonomy of the dis- 
puted species as was recently applicable, for instance, to 
the morphologically uniform bufonid genus Bufotes (Du- 
FRESNES et al. 2019). 

Morphometric indices had a low-resolution power 
for identify hybrids among the Balkan water frog spe- 
cies. For instance, a hybrid MIK3160 possessing 69.9% of 
P. shqipericus genome and 30.1% P. kurtmuelleri clearly fell 
to the centroid of P. kurtmuelleri in a morphological DAPC 
analysis (Fig. 3A.). However, the same individual showed 
an intermediate position in a DAPC analysis of microsat- 
ellite markers, corroborating results of the Bayesian analy- 
sis (Fig. 3B). Similarly, a hybrid MIK3277 possessing two 
thirds of its genome from P. kurtmuelleri and the rest from 
P. epeiroticus was assigned to the P. kurtmuelleri centroid 
in both DAPC analyses. Two other Pelophylax kurtmuel- 
leri x P. epeiroticus hybrids (MIK3233 and MIK3309) also 
showed an intermediate position in a multivariate space in 
DAPC analyses, however, only microsatellite markers sep- 
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arated them clearly from putative pure parental individu- 
als. Therefore, a hybrid identification based solely on mor- 
phometry may lead to an inaccurate estimation of popu- 
lation structure, which is also evident in other amphibian 
taxa (e.g. LEHTINEN et al. 2016, ARNTZEN et al. 2018). For 
detection of hybrids we advocate using nuclear molecular 
markers, including microsatellites and SAI-1. However, the 
length variability of the SAI-1 fragments examined on the 
agarose gel alone may not be sufficient, as shown by the 
example of the P. kurtmuelleri female from Pocgem (Alba- 
nia). This individual possessed SAI-1 fragments of ~720 bp 
and ~840 bp. The first fragment was specific for P. kurt- 
muelleri and the second one resembled the fragment spe- 
cific for P. epeiroticus, indicating that this female is a hybrid 
between these two species. The sampling locality Pocem is, 
however, far northward from the P epeiroticus range limit. 
Sequencing of the ~840 bp fragment and its comparison 
with other sequences using BLAST and a parsimony net- 
work showed that this fragment revealed the highest nucle- 
otide similarity with SAI-1 sequences from P. ridibundus. 
The occurrence of the P. ridibundus-specific SAI-1 allele in 
the genome of P. kurtmuelleri might be a consequence of 
the ancestral polymorphism or historical hybridization be- 
tween P kurtmuelleri and P. ridibundus. Similarly, P. kurt- 
muelleri-specific alleles were found for instance in Poland, 
far away from the supposed range of P. kurtmuelleri (Ko- 
LENDA et al. 2017). 


Conclusion 


The accurate identification of water frogs in the southwest- 
ern Balkans might be a challenging objective due to sub- 
stantial morphological variability, interspecific hybridiza- 
tion, and sympatric occurrence. Despite these challenges, 
the most useful characteristics for the identification of the 
studied Balkan species are morphometric indices L/CINT, 
T/CINT, L/T, and DP/CINT, and qualitative traits like col- 
ouration and the size and shape of the metatarsal tubercle. 
Because some of these field markers may overlap between 
species when used for a taxon identification independently, 
both morphometric and morphological markers are nec- 
essary and recommended. Morphological markers are not 
suitable for the identification of hybrids where the applica- 
tion of genetic markers is necessary. 
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